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! Abstract 

O 

Scattering observables can be computed in lattice field theory by measuring the volume dependence of 

o ' 

r*** ■ energy levels of two particle states. The dominant volume dependence, proportional to inverse powers 



o 



X 



of the volume, is determined by the phase shifts. This universal relation (Liischer's formula) between 
energy levels and phase shifts is distorted by corrections which, in the large volume limit, are exponentially 
suppressed. They may be sizable, however, for the volumes used in practice and they set a limit on how 
small the lattice can be in these studies. We estimate these corrections, mostly in the case of two nucleons. 
Qualitatively, we find that the exponentially suppressed corrections are proportional to the square of the 
potential (or to terms suppressed in the chiral expansion) and the effect due to pions going "around the 
world" vanishes. Quantitatively, the size of the lattice should be greater than « (5fm) 3 in order to keep 
finite volume corrections to the phase less than 1° for realistic pion mass. 
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I. INTRODUCTION 



Recent attempts at studying nuclear interactions using lattice QCD raise an obvious question: 
what is the minimum size of the lattice that can be used in order to accommodate two nucleons 
inside it without significant distortion ? This paper aims at answering this question. 

Lattice field theory calculations are performed using imaginary time. This precludes the cal- 
culation of scattering amplitudes in the infinite volume limit yj]. The usual way of obtaining 
information on scattering amplitudes with lattice techniques is to use the volume dependence of 
the two-particle energy levels yfl. It is intuitively clear that the energy levels of a two-particle 
system are moved up (down) for a repulsive (attractive) interaction and that this shift, due to the 
interactions, vanishes in the infinite volume limit. In general the volume dependence of the energy 
levels is a complicated result of the dynamics but if the size of the box is much larger than the 
range of the interaction between the particles the problem simplifies. In this regime we can sepa- 
rate the volume dependences in two categories: the power law dependences, proportional to 1/L 3 
(where L is the size of the box) and exponentially suppressed contributions proportional to e~ L l R 
(R is the range of the interaction). The power law dependence is fully determined by the elastic 
phase shift y|, |4j, |50 at that energy, through the "Liischer's formula". This relation is universal in 
the sense that it does not depend on the underlying forces between the two particles, only on the 
phase shift at one particular value of the energy. The exponentially suppressed corrections (ESC) 
are less universal. For not too small boxes, however, they are dominated by the lightest particle 
that can be exchanged between the particles; in the case of QCD, those are pions. As those pions 
are soft, standard effective field theory (chiral perturbation theory) techniques can be used to com- 
pute the ESC l . This was indeed recently done, at one loop level, for the case of two pions in the 
i so spin 1 = 2 channel flg]. 

At large enough box volume L 3 the power law dependence is dominant and the exponentially 
suppressed ones can be neglected. In practice however, they may still be sizable and spoil the 
lattice extractions of phase shifts using Liischer's formula. An estimate of these effects in the case 
of two nucleons is particularly urgent now as the first unquenched calculations of nucleon-nucleon 



1 Notice the different level of universality between the power law and the exponential L dependence. The power 
law dependence is the same for two different underlying theories, as long as they have the same phase shift at that 
energy level. The ESC will be the same in two different theories only if they both have the same pattern of chiral 
symmetry and the same low energy parameters. 



2 



FIG. 1: Diagrammatic representation of the Lippmann-Schwinger equation. 

phase- shifts have appeared [|70. In this paper we will discuss how to compute the ESC in general 
and evaluate them numerically using chiral nuclear effective theory in order to estimate their sizes 
in feasible lattice calculations. 



H. LUSCHER'S FORMULA AND EXPONENTIALLY SUPPRESSED CORRECTIONS 
A. Infinite volume scattering matrix 

We first review the relations between the potential, scattering matrix and phase shifts. The 
scattering matrix in the non-relativistic case is given by the sum of the (infinite series of) diagrams 
shown in the first line of Fig. [fl which can be summed up by solving the integral equation shown 
in the second line of Fig. \\\ This integral equation is simply the Lippmann-Schwinger equation 



T(p, k) = -MV(p, k)-M 



d 3 q V(p,q)T(q,k) 



(1) 



(2tt) 3 q 2 - k 2 - iO ' 

where k = \k\ and p = \p\ are incoming and outgoing momenta, and M is the rest mass of nucleon. 
The on-shell amplitude is given by T(k, k) 2 . Projecting onto partial waves we have 

dqq 2 V lmLM (p, q)T LM i>m>(q, k) 



Tlml'm'iP, k) = -MVi m it m >(p, k)- 



LM 



27T 2 



q 2 — k 2 — iO 



(2) 



From now on we'll take the spherically symmetric central potential, i.e., Vimv m ' ~ SwS mm /. In 
the infinite volume limit the partial waves uncouple and from now on we will drop the angular 
momentum indices. 



Our normalization of T includes an extra factor of M compared to the more standard one. 
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The iO prescription makes the amplitude complex and it is important only for q 2 m ME, the 
region of phase space describing the on- shell propagation of the intermediate states over large 
distances. It is intuitively clear that these are the contributions that will receive the largest finite 
volume corrections, as we will see below. It is convenient then to separate from the T-matrix the 
part corresponding to on-shell propagation of intermediate states. For that we rewrite Eq. © as 

TV h\ i/T// ,x vT dqq 2 MV(p,q)T(q,k) . [°° dqq\ 2 2 

T(p,k) = - MV(p,k) - 1 — ^— — i it I —MV{p,q)T(q,k)5{q-k) 

- — - v 

MV(p,k)T(k,k)±- 

.-M,,M,(H#,M,)-ff^ 

where V J represents a principal value integral. Since the inhomogeneous term is multiplied by a 
(p independent) factor, the solution of Eq. ([3]) will be given by 



T(p,k) = -K(p,k) l + i-£-T(k,k) 



(4) 



where K(p, k) satisfies 



The i^-matrix is real below particle production thresholds since its defining equation is real. Notice 
that in the definition of K(p, k) the on-shell propagation of intermediate states is not included on 
the account of the principal value prescription. The on-shell amplitude can now be written as 



T(k,k) = — R -■ (6) 

K(k,k) L 

Comparing with the standard parameterization of T(k, k) in terms of phase shifts S(k) 

Att e i25( - h ^ - 1 47r 

T(k, k) = — - - = (7) 

k 2i k cot S(k) — ik ' 

we find 

We can also write the principal value integral as 

dqq 2 V(p, q)K(q, k) f°° dq q 2 V(p, q)K(q, k) - k 2 V(p, k)K(k } k) 



v 



2tt 2 q 2 -k 2 J 2tt 2 q 2 -k 2 

which is sometimes more convenient, especially in numerical calculations. 



(9) 
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B. Finite volume case 



There are a few changes in the derivation of the Lippmann-Schwinger equation as we go from 
infinite space to a finite box 3 . First, the allowed momenta of intermediate particles are restricted to 
discrete values q = 2nn/ L,n 6 Z 3 . Second, the projection onto partial waves is more complicated 
as the shape of the box breaks rotational symmetry and higher partial waves mix even when the 
potential is spin independent and spherically symmetric. The mixing is however suppressed by 
two effects. The first is that the higher partial waves play little role at low energy scatterings, i.e., 
5i(k) ~ k 2l+1 at small k. In the case of s-wave, for instance, the largest contamination comes 
from the I = 4 partial wave, whose phase shifts are no greater than 2° for < k < 300 MeV. The 
second effect suppressing the partial wave mixing is the approximate orthogonality between the 
spherical harmonics at finite but large volumes. For these reasons, the partial wave mixing will be 
disregarded from now on. Finally, the masses and potentials are also changed from their infinite 
volume values. Those changes are exponentially suppressed and, as we will argue below, are also 
suppressed in the chiral expansion, a point that will be further discussed below. In this section, we 
will keep the infinite volume values of M and V. 

We arrive then at the equation defining T, the finite volume analogue of the T-matrix, 

q=2irn/ L 

Assuming that the value of k does not coincide with any of the allowed values of q there is no 
need for the iO prescription, which makes the finite volume scattering matrix real. We now use 
a fundamental result of finite volume momentum sums [8Q: if the summand has no singularities 
in the real axis and decays to zero fast enough at infinity, the difference between the sum and the 
integral is exponentially small at large L. This result follows from the Poisson summation formula 

q=2nn/L J V ' n^0,n&I? J V ' 

' /,g r/(g)+0(e- mi ), (11) 



(2tt) 3 - 

where m is the characteristic scale of /. The summand in Eq. (flOl) has however a singularity at 
q = k and the Poisson summation formula cannot be directly used. We isolate this singularity by 



3 Of course, there is no real scattering and/or asymptotic states at finite volume, but we will still call the finite volume 
version of T the finite volume scattering matrix. 
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performing similar steps taken above with the infinite volume scattering matrix. We separate the 
singular term corresponding to on- shell propagation by writing 

T( P , k) = -MV(p, *)-■! £ jgkggfa, *) " MV (P- k ™ *) 



q=2irn/ L 



q 2 — k 2 

-MV(p,k)T(k,k)j- 3 J2 



L 3 ' q 2 — k 2 

q=2-Kn/ L 

The first summand is now regular and the sum can be replaced by integrals 

V(p,q)T(q,k)-V(p,k)T(k,k) 



(12) 



M x ^ 
L 3 2s 



q=2irn/ L 



q 2 - k 2 



M 



v f d 3 q V{p,q)T(q,k) mV [ d 3 q V(p, k)T(k, k) 



(2tt) 3 q 2 - k 2 



+ F{p,k), 



(13) 



(2tt) 3 q 2 - k 2 

with the difference between the sum and the integral F(p, k) being an exponentially small quantity 

d 3 q \ V(p, q)K(q, k) - V(p, k)K(k, k) 



k 2 - q 2 



(14) 



The integrand in F is regular, thus the Poisson summation formula can be applied to Eq. (|T4|) to 
yield 

dqq sin(|%L) V(p, q)T(q, k) - V(p, k)T(k, k) 



F(p,k) = -M 



2vr 2 \n\L 



q 2 — k 2 



(15) 



The second sum in Eq. (U2I) includes the power law corrections. The sum itself is a universal 
function of k, independent of the potential. We give it a name by defining 

d 3 q \ 1 



1 s^ i2 



An 2 L V An 2 
We can now write Eq. (fT2l) as 



- y 



q=2irn/ L 



T(p, k) = -MV(p, k) [l + T(k, k) JL_s(^ 



v 



(2n) 3 J q 2 ~ k 2 ' 



d 3 q MV{p,q)T(q,k) 



(16) 



(27T) 



q 2 — k 2 



The solution of Eq. (fTTT) is given by 



T{p,k) = - 1 + T(fc,fc) 



1 „/fc 2 L 2 



Att 2 L \ An 2 

where K(p, k) and ¥(p, k) are, in their turn, defined by 

MV(p,g)% fc) 



K(p, fc), 



F(p, fc). 
(17) 

(18) 



K(p, k) = MV(p,k) - 



v 



q 2 -k 2 



¥{p,k), 



(19) 
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and 

¥(p,k) = -M y r dqq sin ^ L ) V b> g) K fa fc ) - V ^ k ^ fc ) . (20) 

n^ne^ 2712 I"' L q2 ~ k2 

K(p, k) is the a finite volume generalization of the infinite volume K(p, k) that includes the ESC, 
but not the power law finite volume corrections. In fact, dropping the exponentially suppressed 
term F(p, k) from Eq. [191 we recover the defining equation for K{jp, k) (Eq. 0). The poles of the 
"on-shell" amplitude T(k, k) 

nk,k) = - ^ A l (2D 

determine the energy levels in the finite box. 

In the large L limit, the finite volume correction F goes to zero (exponentially fast) and K is 
equal to the infinite volume /^-matrix defined in Eq. ©. Then the relation between the energy 
level shifts and scattering parameters, the Liischer's formula, is recovered: 

4tt 1 k 2 T 2 

-J^k) = kcotS ^ = lTL S ^ <22) 

This formula includes the power law correction to the energy levels in a box, but not the ESC. This 
can be seen more explicitly by approximating, for small k, k cot 5(k) by — l/a (where a is the 
scattering length) and «S(^r) ~ -l/{L 3 k 2 ) to find 



E = — 



k 2 47ra 



(23) 



M ML 3 

If we include F, the difference between the infinite and finite volume matrices, AK = K — K, 
provides corrections to this relation as follows 



k cot 8(k) 



AK(k,k) /AK(k, k) 



+ 



1 „/k 2 L 2 



.L S ll^) (24) 



K(k,k) V K(k,k) 2 
and AK satisfies the equation, 

Jo ( 27r ) Q - k 

Equation (|24l) is one form of the Liischer's formula with ESC associated with the scattering matrix 
due to the finite volume. Further approximations are possible if the pole is close to free particle 
levels, i.e., k 2 ps (2im/L) 2 . 

Let us recapitulate we have done so far. The large (power law) difference between the finite 
and infinite volume scattering amplitude comes from the singular q 2 = k 2 region of the integral. 
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Physically, it corresponds to the kinematics where both intermediate particles are on-shell and 
can propagate far away, explore the lattice and "notice" it is finite. The K-matrix is defined by 
subtracting this term, so the difference from the infinite volume K-matrix is suppressed by terms 
of order e~ mL . The procedure we followed was to isolate the singular term in the finite volume 
T matrix, relate the rest to the K-matrix and the phase shifts and treat separately the singular 
region. The contribution from the singular region, where the power law L dependence resides, 
does not depend on the particular interaction, it is an universal function describing the phase space 
for the two intermediate particles to be on-shell simultaneously and can be computed (numerically 
in general, but analytic formulae are available in limiting cases). This way a relation between the 
phase shifts and the energy levels in a box are found. 

Besides the ESC appearing in the scattering matrix, there are ESC to the potential and particle 
masses themselves. They will be discussed in the next section. 

C. No order V (pions "around the world") effect 

When the nuclear potential V is weak, Eqs. (Q3 and (l20l) show that A¥(k, k) is quadratic in 
the potential. It might seem surprising that there is no ESC linear in V coming from pions being 
exchanged after "wrapping around the world". This can be understood in simple terms. 

In momentum space, the potential generated at finite volume by, for instance, one pion ex- 
change, is the same as in the infinite volume limit. In fact, it is essentially given by the solution of 
the Klein-Gordon equation in momentum space. The potentials in position space however depend 
on the volume. In infinite space we have 

v{rl = I wr^ m - (26) 

but at finite volume we have instead 

q=2ixn/ L 

nezs J 1 ; 

= ^V(r + nL). (27) 

nGZ 3 

The first (n = 0) term in Eq. (ITTT) reproduces the infinite volume result, the remaining ones 
are finite volume exponentially suppressed corrections to it. We see then that there are indeed 
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contributions to the potential between two particles in a periodic box coming from "pions wrapping 
around the world". Let us now compute, in first order perturbation theory, the energy shift due to 
the interaction. For simplicity, let us assume the unperturbed states correspond to plane waves 
with zero momentum. The energy shift is given by 



f 1 1 

AE = J ^ dr x dr 2 — — \{f x - r 2 ) 

=7^ / drV(r) 

J L 3 

JJ / dfV(f+nL) 



l r 

-— / dfV(r) 



h7,,n (28) 



L 3 ' 

where a is the scattering length (at infinite volume). The last line above shows the leading power 
law dependence with the volume contained in the Liischer's formula. We see then that the ESC to 
the potential compensate for the fact that the integrations over space are limited to the finite box 
and the final result, when expressed in terms of the infinite volume quantity a , contains no ESC. 



HI. 1 S NUCLEON-NUCLEON SCATTERING 

In this section we will estimate numerically the size of the ESC in the case of spin-singlet 
nucleon-nucleon interactions. It is important to estimate them not only in the case of realistic pion 
masses but also for the current lattice calculations with pion masses in the = 300 — 500 MeV 
range. This can only be done if the behavior of nuclear forces at large quark masses is known. In 



150. In practice, 



principle, the nuclear chiral effective theory provides this extrapolation [|13L | 
these extrapolations are somewhat hindered by the poor knowledge of some low energy constants 
and the slow convergence of the chiral expansion so the results we present for non-realistic pion 
masses should be taken with a grain of salt. Further lattice QCD/effective theory work in the near 
future should improve the situation markedly. 

In the chiral expansion, the different contributions to the potential can be divided in short 
range (with range on QCD scales r ~ 1/ Aqcd) and long distance (with range on the pion scale 
r ~ l/m^l^l 10, II]. They generate ESC suppressed by factors of, respectively, e~ AQCE,L and 
g-m^L Clearly, the largest one comes from the long distance potential. At leading order in the 
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chiral expansion, the long distance potential is given by one-pion exchange which does not re- 
ceive finite volume corrections since it is defined by tree diagrams. At higher orders in the chiral 
expansion there are contributions coming from two-pion exchange. They decay at large distances 
as e~ 2m "« r and, consequently, generate finite volume corrections proportional to e ~ 2m7lL . The two- 
pion exchange contribution to the energy shifts is thus suppressed by both a factor of e - 2m ^ L and 
a factor of m^/Airf 2 coming from the chiral expansion. Also, at the next order of the chiral ex- 
pansion there are vertex corrections that renormalize the strength of the pion-nucleon coupling 
and that receive (exponentially small) finite volume corrections. We will not include any of these 
higher orders effects to the nuclear forces in our estimates as are they are suppressed in the chiral 
expansion. 

The short distance part of the potential is described, in the effective theory approach, by contact 
terms. The contact terms containing no derivatives form a power series on the quark masses or, 
equivalently, on m 2 : 

C%**' = (C ° + m\Cl + ■ ■ • )(N^t 2 N*)(N t t 2 N), (29) 

where N is the nucleon field and the matrices r act on isospin space. The natural sizes for the 
contact terms are 

°° ~ MA' 

where A is the high momentum scale where the effective theory breaks down, A « 500 MeV. 
Unfortunately, from the experimental values of the phase shift we have access only to the combi- 
nation C = Cq + w^Cq, not to Cq and Cq individually 4 . As such we can only make reasonable 
assumptions based on naive dimensional analysis (Eq. (|30l) ) about the values of Cq and Cq sepa- 
rately, consistent with the value of C at the physical pion mass determined by fitting experimental 
phase shifts. 

The next-to-leading corrections to the short distance potential come from the term 

£ Wo der. = C 2 (N^T 2 V 2 N*)(N T T 2 N) + h.C. (31) 

The long range part of the potential also receives corrections contributing at the same order as 



Processes involving emission/absorption of a pion do distinguish these operators. 
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G%. Both Co and C2 contact terms are fit to reproduce physical phase shifts of nucleon-nucleon 
scattering as is discussed in following subsections. 



A. One-pion exchange potential only 



In order to gain insight and verify our numerics we first consider a case where the ESC can be 
computed analytically. This is the case of a potential including only the one-pion exchange piece. 
Without the short range part the nuclear potential becomes rather weak in this channel and some 
approximations are possible. The 1 S -projected OPEP is given by 

~ml + (q + fc) 2n 



V opc (k,q) 



9a™1 



In 



ml + (q — kf 



(32) 



8f*qk 

where f n = 132 MeV and g& = 1.26. Suppose the potential is weak enough to approximate 
K ~ K ~ MV ope . Then AK ~ F ~ (V opc ) 2 from Eq. (Q3- So, AK/K is roughly proportional 
to V opc . Since | V^ opc | is monotonically decreasing with k, the ESC is largest when k — > 0, which 
is our point of interest. Below we show an analytical form of the ESC associated with the OPEP 
in low energy scatterings. 

We rewrite the F-matrix defined in Eq. (|20l ) in a different form suitable for the following dis- 
cussion of numerical estimates. 



F(p, k) = FH 0°' k )= c ^{P, k), A m ^ 00, 

|m|^0,meZ 3 



n=l 



where 



F m (jp, k) = -M 



dqq sin(|m|Lg) V(p, q)K.(q, k) — V(p, k)K.(k, k) 



(33) 



(34) 



2tt 2 \fh\L q 2 — k 2 

and c n is the number of distinct 3D integer vectors that share a common norm */n. The cubic 
multiplicity c n is provided in Table U for n — 1, • • • ,10. 
The OPEP is expanded in terms of k as 

1 



V ope (k,q) 



2/2 ml + q : 



+ 0(k< 



(35) 



TABLE I: First ten numbers of the cubic multiplicity c n , introduced in Eq. d33l >. 



12345 678 9 10 
6 12 8 6 24 24 12 30 24 
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From Eq. (|20l) , the ESC term F^ corresponding to the lowest order of the OPEP in Eq. (1351 ) is 

M 2 f°° . qsm(^/ELq) 



g\<M 2 



dq- 



dq 



q 2 — k 2 
sin(y/nLq) 
Q 



[(V^(k,q)) 2 -(V°^(k,k)) 2 } 
1 1 



(ml + q 



2\2 



mi 



where in the first line K ~ MV ope approximation is used. The integral can be performed analyti- 
cally and the result is 



dq 



sm(^/nLq) 



(m 2 + q 



2\2 



1 



ml 



71 



2mt 



71 



from formula [3.735] in Ref. 



.HA 



4m* 

Therefore, one obtains 



1 - - e -^ L (2 + Vnm,L) 



71 



2mt 



g\m n M 2 e 



/nm w L 



-(2 + v^m,L) + 0(k 2 



(36) 



327r/* s/nm^L 

Using only the OPEP, we numerically compute AK(k, k) /K(k, k) from Eqs. ©, (fl9i (|20l) . (|25l) 
using the F = 6F 1 approximation, and plot it as a function of L with fixed k <C m^ as shown in 
Fig. [2] In the figure we plot the function R n (L), which has the definition 



n'<n 



/UL)^-^ ' "' 



=MV°p c 



n'=l 
n'<n 



MV°P e (0,0) 



\- c n/ g 2 A m n Me-^ m * L ^ /f 



n'=l 



16nf 2 



n'm-jrL 



(37) 



with n = 1. For a weak potential approximation, AK ~ F and K ~ MV. Therefore the 
function Ri(L) represents an approximated curve for the AK[k, k)/K(k, k) for weak potential 
and low-energy scattering. The exponential suppression e~ m7rL comes from the first term in the 
integer vector summation, X^m^o sm(\fh\Lq), which appears in the definition of the F-matrix. The 
second, third, and n-th terms correspond to corrections proportional to e -V2m^L^ Q-V^m^L^ an( j 
e -Vnrn v L^ reS p ec tively. Finite volume corrections AK(k, k)/K{k, k) using the "partial sums" of 
the F-matrix, J2n=i " CnF^ for finite A n , are plotted in Fig. [3] for A n = 1, 2, 3, 4. Each figure also 
contains a function Rx n (L)- As is clear in the figure, the effects of F^ with higher n are smaller 
than those of the lowest few partial-F, though the convergence of the sum of series is slow. 

The unrealistic case discussed in this section shows that, as expected, the leading ESC are 
proportional to e~ m ^ L . It also reinforces our confidence in the numerics involved. 
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FIG. 2: Plots of AK(k, k)/K(k, k) as a function of L using the OPEP (no contact potentials). Two val- 
ues of gA constant are used: (A) qa = 1.26/vl0 and (B) qa = 1.26. In the case (A) \MV(k,k) — 
K(k, k)\/\K(k, k)\ ~ 0.023, and in the case (B) \MV(k, k) - K(k, k)\/\K(k, k)\ ~ 0.24. In the calcula- 
tion of AK only the F 1 is used as it gives the largest ESC. The cutoff in the potential is set to 500 MeV, the 
pion mass is set to 138.0 MeV, and external momentum is set k = 0.1 MeV The function R\{L) computed 
using the same parameters as the data sets is drawn in each figure. 

B. Realistic potential and pion mass 

In order to estimate the ESC in the realistic case we use the potential discussed above with C 
and C*2 being of order An /MA and An /MA 3 , respectively. There is some latitude in the fitting 
procedure used to determine the short distance constants C and C 2 . In order to explore the 
sensitivity of our results to the choice of constants we selected many sets of (Co,C2), which 
reproduce physical phase shifts for A; up to a few hundred MeV. We determined the value of k 
solving the Liischer's formula (Eq.l22l) and use this value to compute ¥(k, k) and the deviation in 
the phase shift A8(k) due to the ESC, where A5(k) is defined as 

47T 



fccot[5(fc) + A5(k)] 



(38) 



K(k,k) 

The result is plotted, as a function of the box size L in Fig. HI As expected the corrections are 
small and decrease with box size, approximately in an exponential fashion. The ESC shown in 
Fig. H] however do not follow naively expected exponential behavior, ~ e~ mwL (with characteristic 
scale being m w ), but instead, numerical results indicate the characteristic scale being ~ 260 MeV. 
Both k and combine to set the scale for the exponential decay of finite volume corrections with 
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L(fm) 

FIG. 3: Plots of AK(k,k)/K(k,k) as a function of L using the OPEP (no contact potentials). Here 
only the partial sum of the F-matrix Yln=i" c n^^ are used with A n = 1, 2, 3, 4. The potential is weak- 
ened by setting gA = 1.26/\/T0. The cutoff in the potential is set to 500 Me V, the pion mass is set to 
138.0 MeV, and the external momentum is set k = 0.1 MeV. Four corresponding curves are drawn in the 
figure: R X (L), R 2 (L), R 3 (L), R 4 (L). 

increasing box size. To verify this statement we repeated the calculation holding k fixed and small 
(k = 0.1 MeV). As Fig. |5] shows, the agreement between the numerical results and the expected 
g-m^L behayioj- i s roughly observed. 

C. Extrapolation to higher pion masses 

In order to estimate the ESC in current lattice calculations we need to compute them for pion 
masses used in these calculations. The nuclear potential is sensitive to the value of the pion mass 
so, as mentioned before, it is difficult to predict what the nuclear potential at higher pion masses 
is. One might expect however that, as the pion mass grows and the range of the nuclear force 
decreases, the ESC would decrease. In fact, if the value of k coming from the solution of the 
Liischer's formula is used, it sets the scale for the exponential decay of the ESC with L and 
there is little difference between the ESC computed with m n = 138 MeV and m n = 250 MeV. 
This is exemplified in the Fig. [6l The extrapolation is done by first choosing C and C 2 , which 
reproduce empirical phase shifts, and then choosing C® and Cq constrained in the range |Cg \rn* < 
(1/15)|Cq| so that these coefficients are within their naive dimensional analysis ranges. Choices 
of (Co, C2) and ((7°, C5 2 ) are randomly made under the constraints. In this calculation, the values 
of M, f, gA are held fixed at physical values. 
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FIG. 4: The upper figure shows the deviation in the phase shift due to the ESC as a function of box size 
L. Realistic pion masses were used. Contact potential parameters Co and C2 are randomly selected but are 
constrained in such a way that they reproduce empirical phase shifts and stay within dimensional analysis 
range. The F-matrix is approximated to be 6F 1 as it gives the largest ESC (see Eq. (I33T)). The ESC were 
calculated at the value of k determined by the Liischer's formula. The lower figure shows values of the 
lowest (real) fcasa function of L, determined from the Liischer's formula. 

IV. CONCLUSION 

We have derived the relation between energy levels and phase shifts keeping terms that are 
exponentially suppressed in the large volume limit but that can be significant for the lattice sizes 
currently in use. These terms indeed appear suppressed by a factor of e~ mnL as long as they are 
computed at energy much smaller than the pion mass. In actual lattice calculations, the value of 
the energy levels measured are a function of the volume used and are usually in the few hundreds 
of MeV range. That changes the suppression factor: it is still exponential but with a different slope 
than the naive estimate ~ e ~ m7ri . We also show that the effect of pions coming from one nucleon, 
"going around the lattice" and interacting with the other nucleon vanishes in linear order in the 
potential if the energy shift is expressed in terms of the (infinite volume) phase shifts. 
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FIG. 5: Deviation in the phase shift due to the ESC as a function of box size L. Realistic pion masses were 
used and the ESC were calculated at fixed value of k = 0.1 MeV. Contact potential parameters Co and C2 
are randomly selected but are constrained in such a way that they reproduce empirical phase shifts and stay 
within dimensional analysis range. The F-matrix is approximated to be 6F 1 as it gives the largest ESC (see 
Eq. (E3l)). 

We performed numerical estimates using an effective field theory inspired potential. Other 
sources of exponential corrections (shifts in the nucleon mass, pion decay constant, changes in 
the one-loop potential) are argued to be small on the basis of the chiral expansion. Our numerical 
estimates indicate that, in a calculation with realistic pion masses, a lattice size of about (5 fm) 3 
is necessary (and sufficient) for these corrections to stay within 1°. Contrary to naive expectations 
an increase of the pion mass does not substantially reduces the minimum box size necessary for a 
given precision level. 
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FIG. 6: The upper figure shows the deviation in the phase shift due to the ESC as a function of box size L. 
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